Phonon-induced dynamic resonance energy transfer 
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In a network of interacting quantum systems achieving fast coherent energy transfer is a notoriously challeng- 
ing task. While quantum systems are susceptible to a wide range of environmental factors, in many physical 
settings their interactions with quantized vibrations, or phonons, of a supporting structure are the most preva- 
lent. This leads to noise and decoherence inside the network, ultimately impacting the energy-transfer process. 
In this work, we introduce a fundamentally new type of coherent energy-transfer mechanism for quantum sys- 
tems, where phonon interactions actually enhance the energy transfer. Here, a shared phonon interacts with the 
systems and dynamically adjusts their resonances, providing remarkable directionality combined with quantum 
speed-up. We call this mechanism phonon-induced dynamic resonance energy transfer and show that it enables 
long-range coherent energy transport even in highly disordered systems. 

PACS numbers: 87.15.hj,05.60.Gg,71.35.-y,63.20.kk 



Quantum systems exhibit many features that are contradic- 
tory to conventional classical systems. One of the most well 
known is the concept of wave-particle duality, where a quan- 
tum system is able to display both wave-like and particle-like 
properties. Such nonclassical behavior enables quantum sys- 
tems to outperform classical systems for a wide range of tasks 
in information processing and communication 1 1 1. In particu- 
lar, in energy transport, quantum coherence enables a system 
to exploit the dynamics of quantum walks |2|, which are sig- 
nificantly faster than those of classical walks and diff'usion. 
However, despite the many advantages off'ered by quantum 
systems, they are extremely challenging to isolate and invari- 
ably interact with their environment. This causes the loss of 
their wave-like properties by decoherence, ultimately destroy- 
ing any quantum advantage. In this work, we show that when 
quantum systems in a network are designed to interact with 
their environment in a shared way LU, they are able to main- 
tain their quantum coherence and even use it for achieving fast 
coherent energy transfer via quantum walk dynamics that are 
highly robust to disorder. 

We consider quantized vibrational modes, or phonons, 
which are carriers of heat widely regarded as forming one 
of the most dominant mechanisms that reduces coherent ef- 
fects in many physical systems. In energy transport, for in- 
stance, electronic excitations of molecules (excitons) move 
downhill in energy through dipole-coupled molecules, while 
their wave-like behavior diminishes rapidly by relaxation into 
a large number of phonon modes IH. Here, the incoherent 
transfer of an exciton has previously been investigated with 
a shared phonon bath where diff'erent molecules are coupled 
to the same phonon modes |5|. However, recent advances 
in experimental techniques have made it possible to observe 
more complex behavior in the quantum regime for engineered 
structures O and biological systems |7 1. For a shared phonon 
bath, coherent transfer of an exciton has been investigated in 
the equilibrium regime, y » /, where the relaxation rate y of 
phonons to an equilibrium state is faster than the dipole cou- 
pling / between molecules responsible for the exciton trans- 
fer [^. Coherent exciton transfer has also been studied for 



the case that y ^ J where the phonons relax from a nonequi- 
librium to an equilibrium state during the exciton transfer |9|. 
The nonequilibrium behavior of phonons has also been con- 
sidered in the internal electron transfer in a molecule 1 10|. In 
contrast to these earlier studies, here we investigate the ex- 
citon transfer through a chain of molecules in the nonequilib- 
rium regime, y <^ J, where the shared phonon is in a nonequi- 
librium state during the exciton transfer. 

Open quantum systems have their internal energy levels 
shifted when they interact with an environment - the so-called 
Lamb shift, which takes place in a wide range of systems [lT|. 
The shift can even become time-dependent if the environment 
fluctuates temporally, e.g. the AC Stark efl'ect, where oscil- 
lating electric fields cause a temporal shift in energy |[T2ll . 
An example relevant to this work is that of a strongly cou- 
pled system of an exciton and a phonon. Here, the transfer of 
the exciton between dipole-interacting molecules causes the 
phonon to be in a nonequilibrium state, which leads to fluctu- 
ations in the energy levels of the molecules supporting the ex- 
citonic transfer. In this Letter, we investigate a time-dependent 
Lamb shift induced by phonon modes that are shared by the 
interacting molecules, causing the energy-level fluctuations of 
the molecules to become highly correlated. We consider the 
non-adiabatic regime where the frequency of the phonons co 
is comparable in magnitude to the dipole coupling / between 
the molecules. We find that the resonance conditions for the 
exciton transfer are dynamically satisfied as a direct conse- 
quence of this time-dependent Lamb shift due to phonon dy- 
namics. This leads to interesting energy transport features: (a) 
Excitonic hopping between molecules is dynamically condi- 
tioned (dynamic resonance), (b) The exciton transfer can be 
biased toward a particular direction (directionality), and (c) 
The exciton moves rapidly by interfering with itself (quan- 
tum walk). We show that this mechanism, which we call 
phonon-induced dynamic resonance energy transfer (DRET), 
enables long-range energy transport combined with quantum 
speed-up, even in the presence of considerable disorder in 
molecular energies. Our work shows that sharing phonons 
between neighboring molecules not only protects quantum 



coherence |[8j |9l, similar in spirit to decoherence-free sub- 
spaces |[T3l , but also enables the harnessing of quantum co- 
herence for long-range energy transport in highly disordered 
systems. Our results highlight some exciting possibilities for 
future hybrid quantum phononic systems, with the features of 
the described mechanism having the potential to open up new 
functionalities for designing quantum transport devices and 
networks. 

Resonance and disorder. — In both classical and quantum 
energy transport systems, resonance plays a vital role. In clas- 
sical transport, with a phonon environment present, resonant 
incoherent energy transfer occurs for a downhill energy gra- 
dient, a mechanism well known as Forster resonance energy 
transfer (FRET) |i4J . Once an exciton hops incoherently from 
a high to a low energy molecule, the energy difference is dis- 
sipated by phonon relaxation, leading to funneling of the ex- 
citon through the molecules. However, in the presence of lo- 
cal energetic traps, due to disorder in the energy levels of the 
molecules, downhill diffusion leads to trapping of the exciton 
at the energetic minima, suppressing long-range energy trans- 
port. In quantum transport, on the other hand, resonant co- 
herent energy transfer from a donor to an acceptor molecule 
occurs in the absence of energy-level mismatches, as shown 
in Fig. 1(a). In the presence of disorder the transfer then be- 
comes non-resonant and again long-range energy transport is 
suppressed. 

Exciton and shared phonon. — We now show that when dif- 
ferent molecules are coupled to the same phonon, the neces- 
sary resonance conditions for quantum transport can be sat- 
isfied by the dynamics of the phonon, even in the presence 
of significant disorder in molecular energies. Consider the 
case where phonon relaxation is absent and the total energy 
of an excitonic polaron system is conserved. For a chain of N 
molecules, the Hamiltonian of the system in the single-exciton 
manifold is 
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FIG. 1: (Color online) A two-molecule system with initial state 
|/> = |1>. In (a), where (Qi - ^2)//i2 = and fi^fi^ 0, reso- 
nant exciton hopping mediated by a dipolar coupling J 12 occurs and 
oscillates in time. In (b)-(d), where colJn - 1, (^i - Q2) /Ju = 2, 
f\/Jn = 1 and /2//12 = 2, a shared phonon in a nonequilibrium 
state leads to dynamic resonance energy transfer. A quantum pic- 
ture of the phonon dynamics in phase space can be represented using 
the Wigner function lil5J (a video is included in the Supplemental 
Material). 



with 1 7) a single-exciton state where molecule j is excited and 
all other molecules are in their ground states, Ej(p, q) is the 
excitonic polaron energy operator and hJjk is the dipolar cou- 
pling responsible for exciton hopping. The operator Ej{p, q) 
is given by 



^^ 1 2 2 

Ej{p, q) = hQ.j + ^ + 2^^ ^^ ~ ^j) ' 



(2) 



where hQ.j is the excited energy of molecule j and cd is the 
frequency of the phonon mode associated with position q and 
momentum p. The equilibrium position dj of the phonon 
mode is assumed to be deformed by the molecular excitation. 
Using the creation and annihilation operators of the phonon, 
h^ and S, the Hamiltonian can be decomposed sls H = H^ -\- 
Hph-\-He-ph- The first term describes the electronic states of the 
molecules 4 = 1.% H^j + ffoj-') \j} {j\ + Z,-^, hJjk \j} {kl 
the second term describes the phonon H^^ = hcoib^b + |), 
and the last term describes the interaction between exciton 
and phonon 4-ph = - Z^.i Mj \j) {j\ ® (b^ + h with fj = 



^Jm<JJ^(2h)~^dj, leading to fluctuations of the molecule energy 
levels. We assume that the total system is initially in the state 
|/) |0)ph at time ^ = 0, where |/) is the initial single-exciton 
state and |0)ph is the ground state of the phonon defined by 
^ph |0)ph = \h(jj |0)ph. The time evolution of the total system 
is Zt 1 1^> ® l^^(0>, with Pk{t) = (MtMkit)} representing the 
population of molecule k at time t. 

Dynamic resonance. — To clarify the principle of DRET, 
we start with a two-molecule system. Here, once a donor 
is excited (molecule 1), the phonon mode has its equilib- 
rium position displaced by Ji, due to a deformation in the 
structural support of the system. Then the phonon is in a 
nonequilibrium state, moving along an equal energy contour 
Ei(p,q) = Ei(0,0) in phase space of momentum p and po- 
sition q, as shown in Figs. 1(b) and (c). Here, Ei(p,q) is 
the energy of the total system including the electronic energy 
of the donor hQ^i and the phonon energy. £^i(0, 0) is the ini- 
tial energy of the total system. During the phonon 's dynam- 
ics, if the equilibrium position of the phonon mode depends 
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FIG. 2: (Color online) Once an exciton is created at molecule 7, 
the equal energy contours Ek(p,q) = Ej(0,0) are determined in 
phase space. If the contours are close enough such that the to- 
tal energy difference is small compared to the dipolar coupling, 
\Ei(p,q) - E2(p,q)\ < fiJn, near-resonant exciton transfer takes 
place. If not, the exciton transfer is suppressed by non-resonance. 
Here we take colJn = 2.4, (Qi - ^2) I J 11 = 2, /1//12 = -0.5 and 
filJn = 1. In (a), where |/> = |1>, an exciton created at molecule 
1 is transferred from a high to a low energy molecule by dynamic 
resonance. In (b), where |/> = |2>, an exciton created at molecule 2 
is localized at the low energy molecule 2 by non-resonance. 



on the state of the molecules due to a different deformation 
in the structural support, i.e. di 9^ J2, the energy level of 
an acceptor (molecule 2) E2(p,q) will oscillate in time |[T6l . 
While the phonon moves along contour Ei(p,q) = £"1(0,0), 
when the energy level of the acceptor becomes equal to that 
of the donor, Ei(p,q) = E2(p,q), there is a high probabil- 
ity for resonant exciton hopping from the donor to the ac- 
ceptor to take place with the conservation of energy and mo- 
mentum of the total system. This happens at an intersec- 
tion (point (ii) in Fig. 1(c)) between different energy con- 
tours Ei(p,q) = £i(0,0) and E2(p,q) = £i(0,0). When the 
acceptor becomes excited, the phonon moves along contour 
Eiip,^) = £i(0, 0). This leads to temporal oscillations in 
the energy level of the donor Ei(p,q) and non-resonance of 
the exciton transfer until the phonon arrives at another inter- 
section (point (iv)). This phonon-induced dynamic resonance 
energy transfer (DRET) takes place even in the presence of 
electronic energy mismatches Qi ^ Q2. In the presence of 
the mismatch, the radii of equal energy contours are differ- 
ent, as shown in Fig. 1(c), which leads to an asymmetry in the 
time evolution of the populations of the molecules, such that 
the exciton stays at the lower energy molecule 2 for longer, as 
shown in Fig. 1(d). Here the energy of the excitonic polaron is 
partially transferred from the exciton to the phonon and back 
again without dissipation in order to allow the transfer to take 
place. This reversible exchange of energy enables efficient 
energy transfer even in rugged electronic energy landscapes 
with energetic traps, as shown in Fig. SI of the Supplemental 
Material. 
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FIG. 3: (Color online) A seven-molecule linear chain with |/) = |1>, 
Jjk = •/ if 7 = ^ + 1 and Jjk = otherwise. The distance be- 
tween neighboring molecules is assumed to be uniform such that 

m = V<^ = (2Li(^ - l)'^fc(O)'/'. In (a), we take Q,- = Q, 
and fj = for all j and k, which leads to a quantum walk in a 
chain of finite size where an exciton is transferred back and forth 
through the molecules with quantum speed-up due to the reflection 
at the boundaries of the chain. In (b), the chain is coupled to a shared 
phonon mode with nonzero fj. Even in the presence of electronic en- 
ergy mismatches, an exciton can be transferred with quantum speed- 
up by dynamic resonance. When molecule 7, located at one of the 
boundaries, has the lowest electronic energy, the exciton is tempo- 
rally trapped at the molecule due to the non-resonance induced by 
the phonon dynamics. The Supplemental Material provides the val- 
ues of the parameters used and the time evolution of the populations 
of the molecules (cf. Fig. S5). 



Directionality. — Remarkably, DRET can also induce direc- 
tionality, similar to a diode in electrical circuitry. In direct 
contrast to FRET, where the energy transfer is always biased 
toward a lower energy molecule, DRET can induce direc- 
tionality in both downhill and uphill directions, such that the 
transfer becomes biased from a high to a low energy molecule 
Qi > Q2, as shown in Fig. 2, or from a low to a high energy 
molecule Q.i < Q2, as shown in Fig. S3. This directionality 
is induced by the resonance of the transition, such that en- 
ergy transfer in a particular direction is enhanced by dynamic 
resonance, while in the opposite direction it is suppressed by 
non-resonance. Directionality can also be induced in a chain 
of molecules, as shown in Fig. S4. 

Quantum walk. — In a large chain, DRET also enables long- 
range energy transport that is otherwise unavailable. In the 
absence of energy mismatches, a coherent exciton moves very 
quickly by interfering with itself similar to a quantum walk. 
Here the root mean square displacement A(0 = ^j{x^) in- 
creases linearly with time t, A(t) oc t, where x is the distance 
that the exciton moves during time tf2].ln the presence of en- 
ergy mismatches, however, quantum coherence leads to coher- 
ent localization of the exciton, similar to Anderson localiza- 
tion 1 14|. A new and rather surprising feature introduced by 
DRET is that when electronic energy mismatches are present 
and a shared phonon is available, both quantum walk and co- 
herent localization can occur concurrently, where an exciton 
moves through a chain with quantum speed-up and is subse- 
quently localized at a particular molecule for a period without 
losing quantum coherence, as shown in Fig. 3. 

Phonon relaxation and thermal noise. — So far only a sin- 
gle shared phonon mode has been considered. As the number 
of shared phonon modes is increased leading to a shared bath, 
phonon relaxation takes place, suppressing the reversibility of 
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FIG. 4: (Color online) A seven-molecule linear chain with /y^ = 7 if 
] - k±\ and J j]^ - otherwise. In each panel, when an exciton is 
created at molecule 1, it is transferred through the molecules, while 
thermal noise of the phonon bath at temperature T destroys quantum 
coherence, leading to a quantum-to-classical transition in the exci- 
ton transfer. Here, as an example, we take environment couplings 
such that quantum transport takes place up to ? ^ 7iJ~^ and classical 
transport occurs subsequently (cf. Fig. S6). (a) In the presence of 
local baths, when all the molecules have the same electronic energy, 
long-range quantum delocalization takes place due to resonance, but 
no classical funneling occurs due to the absence of a downhill energy 
gradient, A^ ^ 100;r) = V<^ ^ (2Li(^ " l)^^)^^^ - 3.6. (b) On 
the other hand, when the energy gradient is comparable to, or larger 
in magnitude than the thermal energy ^^T, downhill classical fun- 
neling takes place, but non-resonance induces coherent localization, 
leading to short-range quantum delocalization. (c) In the presence of 
a shared phonon bath, however, a time-dependent Lamb shift takes 
place such that a near-flat energy landscape evolves into a downhill 
landscape via phonon relaxation. In this case, long-range quantum 
delocalization and subsequent classical funneling can be utilized for 
enhancing the overall exciton transfer: in (c), the dotted and dashed 
lines display A(0 in (a) and (b) for comparison. The Supplemental 
Material provides a non-Markovian quantum master equation and the 
values of the parameters used. 



DRET. Here, once an exciton hops from a high to low energy 
molecule by dynamic resonance, it is trapped at the lower en- 
ergy molecule due to phonon relaxation, enhancing the direc- 
tionality of the exciton transfer. On the other hand, in the pres- 
ence of thermal noise, DRET is degraded as quantum coherent 
effects are suppressed by decoherence and the trapped exciton 
at the lower energy molecule may be released by absorbing 
thermal energy. Despite the detrimental effects of phonon re- 
laxation and thermal noise, DRET can enhance exciton trans- 
fer in disordered systems by exploiting an interplay between 
quantum delocalization and classical funneling, as shown in 
Fig. 4. This suggests that sharing phonons between molecules 
not only enhances the coherent timescale of exciton transfer, 
as shown in previous studies (IIIH, but also enables the uti- 
lizing of quantum coherence for efficient energy transport in 
disordered systems. 

Summary. — We investigated a shared phonon in the 
nonequilibrium regime and its influence on the exciton trans- 
fer in a chain of molecules. We showed that the time- 
dependent Lamb shift induced by phonon dynamics leads to 
interesting energy transport features, such as dynamic reso- 
nance, directionality and quantum walk. This work opens up 
some intriguing possibilities for quantum transport in disor- 
dered systems in general. Contrary to the widely accepted 
notion that quantum coherence is not useful for energy trans- 
port in disordered systems, our results show that phonons can 
enable a quantum speed-up in energy transport, regardless 



of significant disorder. The consideration of more complex 
shared phonon dynamics may lead to new transport mecha- 
nisms like DRET being discovered, providing features that 
are unavailable in systems investigated so far. Thus, we ex- 
pect our findings to unlock the engineering of novel quantum 
transport schemes in disordered nanostructures. 
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SUPPLEMENTARY INFORMATION 

In section I, we provide the theoretical model of the shared phonon bath used in Figs. 4 and S6. In section II, we provide 
values for the parameters used in all simulations. 

SECTION I. THEORETICAL MODEL OF SHARED PHONON BATH 

Hamiltonian 

We consider a molecular chain coupled to a shared phonon bath modeled by a set of quantum harmonic oscillators. The 
Hamiltonian in the single-exciton manifold consists of three parts H = He -\- H^h + ^e-ph- The first term describes the electronic 
states of the molecules 

where Aj = 2^ ^^/r^J^ ^/l + '^^j) is the Lamb shift induced by the interaction between molecule j and the phonons. The second 
term describes the phonons ^ph = 2^ Hco^^b^, where P^ and b^ are creation and annihilation operators of a phonon mode ^ with 
frequency oj^. The last term describes the interaction between the molecules and phonons 

A^ A^ 

j=l ^ k^l 

When Sj = for all 7, He-ph describes local phonon baths: for a phonon mode ^ coupled to molecule j by gj^ i^ 0, we assume 
g^ = for all k ^ j. For nonzero Sj, H^-ph describes a shared phonon bath where each phonon mode ^ is coupled to all the 
molecules by gj^Sj, except for molecule j which is coupled by gj^(l + Sj). The diff'erence in gj^ leads to decoherence of the 
electronic coherence, while gj^Sj induces a time-dependent Lamb shift without additional decoherence, as shown below. 

The dynamics of the exciton transfer are described by the reduced electronic state pe(0 = Trph[p(0], where p(0 is the total state 
and Trph is the partial trace over phonons. We take the state of the total system at time t = Oto bep(O) = Pe(O)0Z"^ exp(-/57?^ph), 
where Pe(0) is the initial single-exciton state and Z"^ exp(-/57?ph) is a thermal state at temperature T = (kBJS)'^, with Z = 
Trph[exp(-yS^ph)]. 

Derivation of time-dependent Lamb shift 

Here we derive an eff'ective Hamiltonian W^(t) that gives exactly the same exciton transfer dynamics, described by the state 
peif), as obtained by the original Hamiltonian H, 



Pe(t) = Trph[exp(-^^Op(0)exp(^^^0] = Trph 



rexp|-^ r^"V¥^')p(0)r^exp|^ f H'^\f)df 



(3) 



where T and T^ denote the chronological and anti-chronological time ordering operators. The eff'ective Hamiltonian consists 
of three parts W^it) = Hf^{t) -h Hp\^ -h W^ ^, where Hf^{t) is an eff'ective Hamiltonian of the electronic states of the molecules 
including time-dependent Lamb shift 

A^ 

Hf{t) = Y, K^j + ^j + ^Sj Y, g%ojf cos(a;^0) l7) <7l + Z ^^J^ l^'> <^l ' (4) 

and W^ ^ is an eff'ective Hamiltonian describing decoherence of the electronic coherence 

A^ 

^f ph = Z Z ^^^^ i^'> <^'i ^ ^^ + ^^)- (^> 

The effective Hamiltonian provides the analytic form of the time-dependent Lamb shift induced by shared phonon dynamics, 
i.e. IfiSj Yj^g\(ji>7^ cos(6l)^0- It also enables one to employ well-developed non-Markovian approaches |lji2j|, which have been 
developed to describe the interaction between molecules and phonons in the form of Eq. ([5]), as shown below. We derive the 



effective Hamiltonian W^{t) and the time-dependent Lamb shift of Hf^{t) with the use of a path integral representation of the 
reduced electronic state pg(0 and an influence functional developed by Feynman and Vernon |3|. 

To implement a path integral representation, we start with the formal expression of the reduced electronic state of molecules 

Peit) 
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Peit) = Trph 

By splitting the time t into M divisions, i.e. t = MAt, the single-exciton state at time t can be expressed as 
Pe(t) = Trph 
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We then substitute an identity operator 1 = ^f^i 1^) (^1 of the single-exciton subspace at the front and back of each time evolution 



operator exp(-iHAt/fi) or exp(iH^At/fl) for a time increment A^ 
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where the superscripts ± and subscripts / of the kf with non-negative integers / g {0, 1, • • • , M} are introduced to distinguish 
identity operators in difl'erent positions. The partial trace over phonon degrees of freedom Trph in Eq. (fTOl) can be carried out 
with the use of an influence functional [Sj . The reduced electronic state is then expressed as 
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where lik'^^k'^, • • • ,k^,k~f^) is the so-called influence functional that is responsible for the influence of the phonon bath on 
the electronic system dynamics. In the continuous limit, i.e. Ar ^ or equivalently M ^ oo, the influence functional I = 
limA^^o I(kQ,kQ, • • • , Z:^, k^) can be decomposed into two parts 



J = Jz)X Jls, 
where the first term Id is responsible for environmental decoherence of the electronic coherence 



(12) 
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and the second term Jls is responsible for dynamic resonance 



Ils = Y\^w\ -i df [6U, k\f)) - sa k-(m 2s J ^ g%oj-^' cos(oj^f) - Aj 
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with Aj = 2^ 8%^f^^ leading to a time-dependent Lamb shift 2fiSj 2^ S%^^^ cos(oj^f) - AA. Note that the subscript / of kf for 
a finite time increment At is replaced with a continuous variable f or t'' such that kf -^ k-(ti) with ti = I At as At -^ 0. Here, 
^(7, k) is the Kronecker delta defined by S(j, k) = I if j = k and ^(7, ^) = otherwise, and a/r) is a bath response function 



characterized by the spectral density Aj(oj) = 2^ ^g%^(p^ - ^^) 
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Hence, in the continuous limit Ar -^ 0, the reduced electronic state can be expressed as 
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with ti = I At, Hf^if) is the effective Hamiltonian of the electronic states of the molecules including time-dependent Lamb shift 
in Eq. ^. This implies that the effective Hamiltonian W^{t) gives exactly the same exciton transfer dynamics peit) as obtained 
by the original total Hamiltonian H. For an arbitrary form of the spectral densities Aj(oj), the dynamics of the reduced electronic 
state pe(t) can be calculated using a path-integral approach, known as the iterative real-time quasiadiabatic propagator path- 
integral (QUAPI) scheme |T| . In this work, however, we take an Ohmic spectral density with a Lorentz-Drude cutoff function 
Ajioj) = Z^ %y/r^(^ - ^^) = (2'hAjln)(jjyjl{(j? -h y^^ with the phonon relaxation rate jj to employ the so-called hierarchical 
approach, which is a computationally efficient framework to simulate exciton transfer dynamics |2|. The time-dependent Lamb 
shift then becomes IfiSj Z^ g^-^<^'^^ cos(^^0 = 2fiSjAje~^J^. This implies that the shared phonon bath induces dynamic resonance 
such that the effective energy landscape that an exciton feels during transport changes in time fi^j -\- fiAj(l -h 2sje~^J^). 

Hierarchically coupled master equations in the presence of a shared phonon bath 

Here we provide a non-Markovian quantum master equation, represented by a collection of hierarchically coupled master 
equations, that reliably describes exciton transfer dynamics in the presence of a shared phonon bath. 

With the shared phonon bath modeled by the Ohmic spectral density and the high temperature condition characterized by 
pfiyj < 1 taken |2|, the time evolution of the reduced electronic state pe{t) can be modeled by the following hierarchically 
coupled master equations 

r . A^ A^ A^ 

j-i j-i M 

where Hf^{t) is the effective Hamiltonian of the electronic states in Eq. ^ and n = {^i, • • • ,^a^} with non-negative integers 
Hj e {0,1,'" }. The vectors nj+ differ from ft only by changing rij to rij ±1, i.e. nj+ = [ni,- • • ,nj ±\,- • • , n^}. In the above, the 
reduced electronic state is the lowest rank operator pg(0 = ^(0, t) with = {0, • • • ,0}. The other auxiliary operators &{n, t) are 
introduced in order to take into account environmental decoherence of electronic coherence. Their dimensions are the same as 
that of the reduced electronic state pe(0. At initial time r = 0, just after an exciton is created in a molecular chain, a-(0, 0) = Pe(0) 
and all the other auxiliary operators are set to be zero, i.e. &(n, 0) = for all ft ^ 0. The last two terms in the coupled master 
equations produce a coupling between the operators of different ranks A^rank = Z^i ^> which are given by 

^jl&iftj^, t)] = i[\j) {j\ , &{ftj^, 0], (19) 

^jV&{ftj., 0] = i^^^ilj) {j\ , ^(ftj-, m + Ajyj{\j} {j\ , &{ftj., t)]. (20) 

The ranks of the auxiliary operators continue to infinity, which is impossible to treat computationally. In order to terminate the 
hierarchically coupled master equations, we set the higher rank operators of A^rank > A/^cutoff to be zero, i.e. &{ft, = for ft 
satisfying A^rank > A^cutoff • In numerical simulations, this means we increase A/^cutoff until the dynamics of the reduced electronic 
state show convergence. In this work, we take A^cutoff = 10 in the simulations of Figs. 4 and S6. 

We note that the hierarchically coupled master equations in Eq. ( 18 ) are different from those in Ref. fT|, where local phonon 



baths were considered. The local bath model corresponds to the case that Sj = for all j where the effective Hamiltonian of the 
electronic states of molecules Hf^(t) is reduced to the original Hamiltonian H^ as the time-dependent Lamb shift vanishes. 

SECTION II. PARAMETERS USED IN SIMULATIONS 

Here we provide values for the parameters used in the simulations. 

Dynamic resonance induced by a shared phonon mode 

In Fig. 1, we consider a two-molecule system with N = 2 and take the initial single-exciton state |/) as being in a pure state 
|1) where molecule 1 is excited and molecule 2 is in its electronic ground state. In Fig. 1(a), we take 1) (Qi - 0.2)1 Jn = (no 



electronic energy mismatch between molecules), and 2) /i = /2 = (no interaction between molecules and phonon mode). In 
Figs. l(b)-(d), we take 1) (Di - ^2)1 Jii = 2 (electronic energy mismatch), 2) co/Ju = 1, fi/Ju = 1 and /2//12 = 2- here the 
dipolar coupling Ju is small in magnitude compared to the electronic energy mismatch, which would lead to non-resonance of 
exciton transfer if the shared phonon mode was not coupled to the molecules. 

In Fig. SI, we consider a three-molecule system with N = 3 and take the initial single-exciton state |/) as being in a pure state 
|1). We take 1) Ju = J23 = J and /13 = (a linear chain), 2) (D2 - ^i)/J = 3 and (Q3 - Qi)// = 1 (a rugged uphill electronic 
energy landscape, Q2 > ^3 > ^1), implying that molecule 2 is an energetic barrier and molecules 1 and 3 are nonresonant. In 
Fig. SI (a), we take /i = /2 = /s = (no interaction between molecules and phonon mode). In Fig. Sl(b), we take oj/J = 1, 
fi/J = 2.11, f2/J = 2.80 and fs/J = 2.56. 

In Figs. 2 and S2, we consider a two-molecule system with N = 2 and take 1) oj/Ju = 2.4, 2) (Qi - 0.2)! Jii = 2, 3) 
fi/Ju = -0.5 and /2//12 = 1. In Figs. 2(a) and S2(a), we take the initial single-exciton state as being in a pure state |/) = |1) 
where molecule 1 is excited and molecule 2 is in its electronic ground state. In Figs. 2(b) and S2(b), on the other hand, we take 
|/) = |2) where molecule 2 is excited instead. 

In Fig. S3, we consider a two-molecule system with N = 2 and take 1) oj/Ju = 0.143, 2) (Qi - ^2)/Ju = -2.14, 3) 
fi/Ji2 = -0.857 and /2//12 = 0.143. In Fig. S3(a), we take the initial single-exciton state |/) = |1), whereas in Fig. S3(b), we 
take |/) = |2). 

In Fig. S4, we consider a three-molecule linear chain, with N = 3, Ju = J23 = J and Ju = 0, and take 1) (Qi - ^3)/ J = 4 
and (D2 - ^3)/^ = 2 (a downhill electronic energy landscape), 2) oj/J = 2.29, fi/J = -0.975, f2/J = 0.654 and fs/J = -0.654. 
In Figs. S4(b)-(d), we take the initial single-exciton state as being in a pure state |/) = |1), |/) = |2) and |/) = |3), respectively. 

In Figs. 3 and S5, we consider a seven-molecule linear chain, with N = 7, J jk = J if j = k±l and Jjk = otherwise, and 
take the initial single-exciton state |/) = |1). In Figs. 3(a) and S5(b), we take 1) (Qy - ^k)/J = for all j and k (a flat electronic 
energy landscape), and 2) fj = for all j (no interaction between molecules and phonon mode). In Figs. 3(b) and S5(d), we take 

1) (^i - Clj)/J = 1.93, (^2 - ^i)IJ = 2.06, (D3 - D7)// = 2.1 1, (Q4 - ^7)// = 2.13, (Q5 - ^7)// = 2.14, (Qg - ^7)// = 2.05, 

2) oj/J = 1, fi/J = -0.471, f2/J = -0.305, fs/J = -0.221, 74// = -0.151, fs/J = 0.129, fe/J = 0.325 and fi/J = 1.47. 

Dynamic resonance induced by a shared phonon bath 

In Figs. 4 and S6, we consider a seven-molecule linear chain with N = l, J jk = J if j = k±l and Jjk = otherwise, and take 
the initial single-exciton state as being in a pure state pe(0) = |1) (1|. In Figs. 4(a) and S6(c), we take 1) (flj - ^k)/J = for all 
j and k (a flat electronic energy landscape), 2) Aj/J = 0.35, yj/J = 0.35, ksT/fiJ = 2 and Sj = for all j (local phonon baths). 
In Figs. 4(b) and S6(d), we take 1) (Q^ - Q^+i)// = 2 (a downhill electronic energy landscape), 2) Aj/J = 0.35, yj/J = 0.35, 
ksT /fiJ = 2 and Sj = for all j (local phonon baths). In Figs. 4(c) and S6(e), we take 1) (flj-flj+i)/J = 2 (a downhill electronic 
energy landscape), 2) Aj/J = 0.35, yj/J = 0.35, UbT /fiJ = 2 and Sj = j x (kBT/HAj) (a shared phonon bath). 
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(a) Rugged electronic energy landscape (no shared phonon) 
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(b) Rugged electronic energy landscape (shared phonon) 
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FIG. S 1 : Energy transfer through a rugged uphill electronic energy landscape mediated by reversible exchange of energy between exciton and 
phonon. (a) In the absence of a shared phonon, an exciton is localized at the initial molecule in a rugged landscape due to the non-resonance 
arising from the energy-level mismatches between initial molecule 1 and the other molecules. Here, a rugged uphill electronic energy landscape 
Q2 > ^3 > Hi is considered, where the energy-level mismatches are similar in magnitude to the dipolar couplings / between the molecules, 
(b) In the presence of a shared phonon, an exciton can be efficiently transferred through the rugged landscape. When the exciton hops from 
low energy molecule 1 to high energy molecule 2, the electronic energy is increased (Qi < ^2), while the potential energy of the phonon is 
decreased \m(j?{q - dif > ^ma/(q - ^2)^ due to the change of the equilibrium position of the phonon mode (di -^ ^2)- In an opposite way, 
when the exciton hops from high energy molecule 2 to low energy molecule 3, the electronic energy is decreased (Q2 > ^3), while the potential 
energy of the phonon is increased. With energy conservation satisfied, quantum coherence then enables a unidirectional exciton transfer from 
molecule 1 to 3, such that the exciton stays at molecule 3 for a while with a high probability. 
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(a) Dynamic resonance energy transfer from high to low energy molecule 
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(b) Non-resonant energy transfer from low to high energy molecule 
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FIG. S 2: Directionality of downhill energy transfer induced by dynamic resonance. Fig. 2 is displayed in detail with the excitonic polaron 
energy surface. 



(a) Dynamic resonance energy transfer from low to high energy molecule 
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(b) Non-resonant energy transfer from high to low energy molecule 
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FIG. S 3: Directionality of uphill energy transfer induced by dynamic resonance. Here we consider a two-molecule system with Qi < Q2 
where electronic energy of molecule 1 is small in magnitude compared to that of molecule 2. In (a), an exciton created at molecule 1 is 
transferred from low to high energy molecule by dynamic resonance. In (b), for the same choice of system parameters, an exciton created at 
molecule 2 is localized at the initial high energy molecule by non-resonance. 
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(a) A three-molecule chain coupled to a shared phonon mode 
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(b) Dynamic resonance energy transfer from molecule 1 to molecule 3 
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(c) Dynamic resonance energy transfer from molecule 2 to molecule 3 (non-resonance between molecules 1 and 2) 
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(d) Non-resonant energy transfer from molecule 3 to molecules 1 and 2 (non-resonance between molecule 3 and the other molecules) 
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FIG. S 4: Directionality of dynamic resonance energy transfer in a chain of molecules. Here we consider a three-molecule system with a 
downhill electronic energy landscape, Qi > Q2 > ^3 as shown in (a). In (b), an exciton created at molecule 1 is resonantly transferred from 
molecule 1 to molecule 3 by dynamic resonance. In (c), an exciton created at molecule 2 is resonantly transferred to molecule 3 by dynamic 
resonance, while exciton transfer to molecule 1 is suppressed by non-resonance. In (d), an exciton created at molecule 3 is trapped at the initial 
molecule with high probability due to the non-resonance between molecule 3 and the other molecules. 
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(a) A linear chain of molecules 
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(c) A linear chain of molecules coupled to a shared phonon mode 
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FIG. S 5: Quantum walk and coherent localization induced by a shared phonon mode. In (a), a linear chain of molecules is shown with 
uniform dipolar couplings / between neighboring molecules. In (b), the dynamics of the populations of the molecules are shown, which leads 
to the root mean square displacement A(t) - ^J{x^) shown in Fig. 3(a). In the absence of electronic energy mismatches, a resonant quantum 
transition takes place between the molecules. This leads to a quantum walk in a chain of finite size where an exciton is transferred back and 
forth through the molecules with quantum speed-up due to the reflection at the boundaries of the chain. In (c), a linear chain of molecules is 
coupled to a shared phonon mode. In (d), the dynamics of the populations of the molecules are shown, which leads to the root mean square 
displacement A{t) = ^/{x^) shown in Fig. 3(b). Even in the presence of electronic energy mismatches, an exciton can be transferred with 
quantum speed-up by dynamic resonance. When molecule 7, located at one of the boundaries, has the lowest electronic energy, the exciton 
is temporally trapped at the molecule due to the non-resonance induced by the phonon dynamics. This is a valuable characteristic to have for 
quantum transport through a chain, for which a source is coupled to one of the boundaries (molecule 1) and a drain (or sink) is coupled to the 
opposite end (molecule 7). When the internal couplings of the chain, such as the dipolar couplings, are large in magnitude compared to the 
external coupling between the chain and the drain, the temporal energy trapping at the boundary will enhance unidirectional energy transport 
from the source to drain by suppressing the reflection at the boundary for a period. This coherence-preserving energy-trapping scheme may 
also be useful for the transfer of phase information, such as in a quantum state. 



(a) A chain of molecules coupled to local phonon baths 
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(b) A chain of molecules coupled to a shared phonon bath 
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(c) Local phonon bath, 
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(e) Shared phonon bath, 
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FIG. S 6: Long-range quantum delocalization and classical funneling induced by a shared phonon bath. In (a), a linear chain of molecules is 
shown, where each molecule is coupled to an independent phonon bath. In (b), a linear chain coupled to a shared phonon bath is shown. In 
(c)-(e), the dynamics of the populations of the molecules are shown, which correspond to the root mean square displacement A(0 = V(^ 
displayed in Figs. 4(a)-(c), respectively. 



